Motion of a condensate in a shaken and vibrating harmonic trap 
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The dynamics of a Bose-Einstein condensate (BEC) in a time-dependent harmonic trapping 
potential is determined for arbitrary variations of the position of the center of the trap and its 
frequencies. The dynamics of the BEC wavepacket is soliton-like. The motion of the center of the 
wavepacket, and the spatially and temporally dependent phase (which affects the coherence prop- 
erties of the BEC) multiplying the soliton-like part of the wavepacket, are analytically determined. 



Bose-Einstein condensates of dilute atomic gases in magnetic traps provide a simple many-body system in which to 
investigate the evolution of a macroscopic coherent quantum system under the influence of external forces . Analytic 
mean- field solutions for these systems exist even for time dependent external forces |2|,|| . The special scaling properties 
of the harmonic potential, created by the interaction of the atomic magnetic moments and the average magnetic field 
of the trap which has a quadratic spatial form, makes it easy to determine the evolution of the condensate even under 
temporal variations of the frequency of the trap f§,||. Moreover, Kohn || and Dobson |?j have shown that, for any 
many-body system in an arbitrarily changing harmonic potential, the motion of the center-of-mass of the system is 
decoupled from the motion of other degrees of freedom of the system. 

Here we present an exact solution for the motion of a Bose-Einstein condensate under the influence of a harmonic 
magnetic field whose center moves as an arbitrary function of time and whose frequency varies arbitrarily with time. 
When the frequency of the harmonic trap is constant in time, the motion of the condensate is as a rigid body whose 
shape is not changed as the potential moves, the motion of the center of the condensate is analytically determined, 
and the time dependence of the phase of the condensate is also analytically obtained. We stress that this result 
applies not only at the mean-field level of approximation, as described by the Gross-Pitaevskii equation [gj , but quite 
generally at the field-theory level. Hence, upon shaking a harmonic potential, no matter how vigorously or quickly, a 
condensate docs not develop an above-the-mean-field component, and does not develop a temperature. No amount 
of shaking will yield a thermal cloud in such a system. When the frequency of the harmonic potential also varies, the 
center of mass motion of the condensate and its phase can still be analytically determined, and the condensate shape 
is obtained by solving the Gross-Pitaevskii equation for a harmonic potential whose center is not moving. 

We consider a general system of N mutually interacting identical Bosonic particles of mass m in an external time- 
dependent harmonic potential. The Hamiltonian for the system is given by 



H = J2 



t—x : y,z 



- Ufa-**). (1) 



Here XQi (t) is the ith component of the time-dependent position vector of the center of the harmonic trap and x a i 
is the ith component of the evth Boson. Let us make the coordinate transformation to a new system of coordinates 
comprising X = ^ a x a /7V and q a = x Q — X for a = 1, N — 1. Note that q^r = x^v — X is a dependent variable 

equal to — ^2a=i Q " an d 12a=i c l a = 0- We can write the quadratic term of the harmonic potential using these 
variables as: 



~m Y ^i{t)[Y^i + Xi-x Qi {t)f= Y ">?(*) [l M {X l -x 0l (t)} 2 + ^mY' 

i—x,y,z a i—x.y,z \ a 



, _ Qai I , (2) 

\ a 

where we have defined M = Nm and we used the relation J^ a=1 q a i — 0. If, for simplicity of notation, we drop the 
indices i that specify the different components of the harmonic frequency and position variables, the above Hamiltonian 
now reads 



JV-l 



H = [P 2 /2M + ^-Muj 2 (t)(X - x {t) 2 } + Y 



^- + l - m u 2 {t) q l + Y U ^-^ 



(3) 



1 



Hence, the general wavefunction tjj(xi, X2, xn ,t) may be written in product form, ip(xi,X2, ■■■,X]y,t) = 
r >l>CM(X,t)ip re i(qi,q2,...,qN-i), where the relative part of the wave function does not depend explicitly on time, 
the effect of the motion (time dependence) of the harmonic potential is only on the center of mass part of the wave 
function, and this dependence is given via the the quantity xo(f). Thus, the motion relative to the center of mass is 
decoupled from the motion of the center of mass, and only the latter is influenced by the shaking. This is true at the 
field-theory level. 

At the mean-field theory level, the harmonic potential appearing in the Gross-Pitaevskii equation is of the form 



V(x,f) = 



i— x,y,z 



Ui{t) 2 [x, - x M (t)} 2 = V&(x-xo(t),t) 



(4) 



where xo(t) = (xo(t), yo(t), Zo(t)) is an arbitrary time-dependent displacement of the center of the trap and 
(uj x (t),u)y(t),uj z (t)) are the (perhaps time-dependent) trap frequencies. Under the action of a potential of the form 
(||) , the internal dynamics of a system of particles is not affected by an arbitrary motion of the center of the potential. 
This follows from the fact that a quadratic potential of the form (||) can be expanded at any instant of time around 
the mean value R(t) = (R x (t), Ry(t), R z {t)) of the center-of-mass of the system as 



V{x,t) = Vb(x-R) 



Xio)Xi + — }_^uj 1 {XiQ-Ri) 



(5) 



The second term in Eq. (|^) corresponds to a homogeneous time-dependent force acting on the system of particles and 
it can therefore be responsible only for a global shift of the wavefunction in position and momentum space. The third 
term is coordinate-independent and can therefore be responsible only for a global phase shift of the wavefunction. 

We demonstrate this general principle by considering a condensate of alkali atoms in a time-dependent harmonic 
magnetic trap as given by Eq. (Q); we produce the general form for the condensate wavefunction. The following is 
thereby a generalization of the result of Heller Q for time-dependent harmonic potentials to the case of interacting 
many-body boson systems. For simplicity we assume that the condensate can be described by a single mean-field 
wavefunction ijj(x,t), but this assumption is not necessary because the following treatment can be automatically 
applied if ip(x,t) is regarded as a field operator. In mean-field, the dynamics of a condensate of weakly interacting 
atoms is determined by the time-dependent Gross-Pitaevskii equation, 



ih 



~dt 



h'v 2 

2m 



+ V(x, t) + U a \il>(x,t)\ 



(6) 



where U = —N ao is proportional to the number of atoms in the condensate and ao is the s-wave scattering length 
for collisions between the atoms. Under the influence of the potential of the form (0), the solution of this equation 
can be written as 



V>(x, t) = Vo(x - R(t), t) exp {*[P(t) • x/ft - 0(t)]} , (7) 

where ^>o( x ,t) satisfies the time-dependent Gross-Pitaevskii equations with Xo(t) = 0, and the ith component of the 
vector R(t) satisfies the equation of motion 

R t + {oj l (t) f(R l -x l0 {t)) -0 . (8) 

In Eq. (Q) the momentum P(t) = mR(t), and the phase <p(t) = J2i 4>i{t) is given by 

^ = ^T,{l^i(t')[RUt')-^o(t')} + [ik{t')fdt^ . (9) 

The proof is straight-forward. Assume that ^(x, t) satisfies Eq. (||) with xo = 0, then substitute the solution (Q) 
into Eq. (^) for arbitrarily varying xo (t) to obtain 

-mR'X + fi^= yR 2 + ^mJ2 UJ ii(^ ~ x mf - {xi - R l ) 2 ]^o . (10) 

i 

It is easy to verify that this is indeed satisfied for Ri(t) and <f>(t) given above. 
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Note that the phase factor exp{z[P(i) • x/7i — 0(t)]} in Eq. (Q) affects the coherence properties of the condensate 
wavepacket, which are given in terms of the coherence function C (p, r;t) = J d 3 x ip* (x + p, t + t)^(x, t) ||,[l0| . 

As an example, let us first consider the solution in the case where u> x (t) — uj x = const and determine R x (t) and 
4> x {t). In this case ipo (x) may be taken as the ground-state solution of the time- independent Gross-Pitaevskii equation, 
which is obtained by replacing the time derivative on the left-hand side of Eq. (||) by ptpo, where p is the chemical 
potential of the condensate. In this case the solution (^) is a solitary solution, namely, the condensate moves as a 
rigid body without changing its shape. The general solution for R x (t) in Eq. (|J) is then given by 

'in a; t 

R x {t) = cosa^tR^O) + bmUJx R x (0)+oj x / dt' smw x (t - t')x (t') . (11) 



In what follows we give some examples for the dynamics in different shaking schemes. For a periodically shaken trap, 
such that xo(t) — xo smut and R x (0) = 0,R X (0) = we obtain the following solution 

Rx(t) = tt-^ o [^xSinujt - lu sin u} x t]x , (12) 

0J X — U) A 

so the instantaneous difference of the center of the wavepacket from the center of the potential is given by 

R x (t) — xo(t) = —5 ^ [uosinuot — u> x sinuj x t] xq . (13) 

LO x — UJ Z 

The expression for the phase 4> x (t) can be easily obtained analytically in terms of simple trigonometric functions using 
Eq. (|). 

We identify three different regimes for this solution: 

1. The adiabatic regime: When u> <C w X) the motion of the condensate will adiabatically follow the motion of 
the center of the trap. We can then approximate 

Rx(t) ~ x o(t) ~ — [ — sin cut — sin uj x t]xo , (14) 

U) x U) x 

whose maximal value is max{\R x (t) — x (t)\} w xqlo/lo x . The adiabatic approximation is justified if this is much 
smaller than the spatial width Ax of the condensate wavefunction V'Oj i-e., lu/lu x <^ Ax/2xq. 



2. Resonance: If u = ui x , then 



R x {t) = [sinuo x t — u x t cosu) x t] , (15) 



and the amplitude grows linearly with time. 
3. Averaged effective potential: When ui ^> u> x , 

R x (t) rs xq— sinuj x t , (16) 

and the motion of the center of the condensate is only slightly affected by the shaken trap. 

Fig. plots R x (t) versus t for four different values of the ratio oj/uj x . In Fig. [pa, oj/u> x = 0.3 and the small deviation 
of R x (t) from xo(t) is evident. In Fig. |l|b, oj/uj x = 0.8 and a substantial overshoot of R x (t) relative to xo(t) is obtained. 
The linear growth of R x (t) is clearly seen for oj/ui x — 0.99 in Fig. [j]c, and the small oscillation of R x {t) is evident for 
lo/oj x = 5 in Fig. [j]d. 

For two dimensional motion corresponding to an ellipsoidal rotation of the center of the trap, xo(t) = xq sin tot and 
2/0 (t) = j/o [ 1 — cos ujt], the solution for R x is given by Eq. ( |l2] ) and R y is 

„ . , 1 — COS LOi,t LO v , , 

RJt) = v — + y n [cos LOyt - cos cut] . 17 

UJy LOy UJ 

The phase cj>(t) is again simple to calculate. 

If the trap is periodically shaken in the x direction for a finite duration to and then the shaking is stopped, then 
the solution for the center of the wavepacket for < t < to is given by Eq. (fi"^), and for t > to, 
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R x (t > to) = — — - — -xq {[lu x sinw^o — ^ sin o> x io] cosui x t + u> [coswio — cosu x to] sinu> x t} . (18) 

Analytic solutions to Eq. (ra) for arbitrary uoi{t) are not known. Even for harmonically varying u>i(t) and Xio{t), 
where Eq. (||) corresponds to a driven Mathieu equation JD]], analytic solutions are not available. In this case, 
depending on the ratio of the frequencies for the variation of u>i(t) and Xio(t), regular bounded motion or unbounded 
motion of Ri(t) may result. Nevertheless, Eq. ([?]) with central wavepacket coordinate R(i) and phases 4>i{t) given by 
Eqs. (Q) @ gives the analytic form for the wavefunction. 

The following numerical example illustrates the nature of the analytic solution for a shaken vibrating trap. Fig. |^ 
plots R x (t) versus t for a periodically varying trap frequency Lo x {t) = o; x (0)[l + 5 sm(uj v t)]. where 5 and lo v are the 
amplitude and frequency of the vibration. Unbounded motion is expected when either the shaking or the vibration 
frequencies are resonant with the basic trap frequency ui x (0) (Fig. |^b), or if the sum of the vibration frequency and 
shaking frequency, or the sum of integers times these frequencies, is resonant with u> x (0) (Fig. j2^,d). Unbounded 
solutions for the volume of the condensate are expected even in the absence of shaking when the vibration frequency 
is resonant with the trap frequency ^,^|. However, the motion of the center of the trap causes the center of mass of 
the condensate to move with respect to the center of the trap, and the amplitude of this motion may be amplified by 
a resonant change of the trapping frequency. 

We stress that the main result of this paper is valid for any system of interacting particles in a harmonic potential. 
Multi-component BECs with harmonic potentials have solutions which also have the properties discussed above, 
provided that the time-dependent harmonic potentials for the various components are exactly the same. Since the 
magnetic moments of atoms in different Zeeman levels, or different isotopic species of the same element, or of different 
elements (for mixed species BECs) are in general not identical, our solution will in general not be relevant for these 
cases. 

In practice, magnetic traps for BEC alkali atoms are harmonic only near their center within a range r/,, which 
is typically a few tens of microns. If such a magnetic trap is shaken, the condensate may enter a region where the 
true potential is anharmonic. In this case the shape of the condensate may change and the motion of its center 
of mass may deviate from the solutions given above. The maximal velocity that can be achieved by accelerating 
a condensate within the harmonic range is roughly given by v max = uJtrapfh, which is typically of the order of 
a few cm/sec. By shaking a magnetic trap using a time-dependent gradient magnetic field it is possible to boost 
condensates in a controlled manner to this range of velocities without changing their shape. This method can also be 
used in conjunction with other methods of optically output coupling high momentum wavepackets p^JT^ | to create 
novel wavepackets. 

To summarize, quite generally, shaking a a harmonic trap will not cause a thermal cloud of atoms to develop from 
a condensate state; only the center of mass motion of the condensate is affected. This is true even at the field theory 
level. We determined analytic solutions for the dynamics of Bose-Einstein condensates of dilute atomic gases in shaken 
and vibrating harmonic traps as described by the Gross-Pitaevskii equation (the mean-field level of approximation). 
One potential application of relevance for the field of atom optics is to boost BECs to desired velocities. 
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No. 212/01), the Israel MOD Research and Technology Unit, and the National Science Foundation through a grant 
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FIG. 1. R x (t) for a periodically shaken trap with constant trap frequency lo x and four different values of the shaking frequency 
u>. The dashed curve is the position of the center of the trap xo(t). 
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